Disorder induced transition between s± and s ++ states in two-band superconductors 
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We have reexamined the problem of disorder in two-band superconductors, and shown within the 
framework of the T-matrix approximation, that the suppression of T c can be described by a single 
parameter depending on the intraband and interband impurity scattering rates. T c is shown to 
be more robust against nonmagnetic impurities than would be predicted in the trivial extension of 
Abrikosov-Gor'kov theory. We find a disorder-induced transition from the s± state to a gapless and 
then to a fully gapped s++ state, controlled by a single parameter - the sign of the average coupling 
constant (A). We argue that this transition has strong implications for experiments. 
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Introduction. The symmetry and structure of the su- 
perconducting order parameter in recently discovered 
iron-based superconductors (FeSC) is one of the main 
challenges in this field^. The Fermi surface (FS) is 
usually given by two small hole pockets around the 
T = (0, 0) point and two electron pockets around the 
M = (jr, 7r) point in the 2-Fc Brillouin zone. The proxim- 
ity of the competing spin-density-wave (SDW) state with 
Q = (71", 7t) suggests antiferromagnetic fluctuations as a 
mechanism for electron pairing. In this case, the natural 
order parameter for most of the FeSC is the so-called s± 
state, described by an isotropic order parameter on each 
FS with the opposite signs for electronlike and holclikc 
pockets. Many experimental results, such as the NMR 
spin-lattice relaxation rate, the spin-resonance peak at 
the SDW wave vector Q in inelastic neutron scattering, 
and quasiparticle interference in tunneling experiments, 
are in good qualitative agrement with this scenario, al- 
though some materials are more anisotropic than others 2 -. 

Since varying amounts of disorder are present in the 
materials, and because superconductivity is created in 
some cases by doping, it is important to understand the 
role of impurities. It has been shown that in an s± state, 
any nonmagnetic impurity which scatters solely between 
the bands with a different sign of the order parameter 
suppresses T c in the same way as a magnetic impurity in 
a single-band BCS superconductor—. Therefore the crit- 
ical temperature T c should obey the Abrikosov-Gor'kov 
(AG) formula In T c0 /T c = #(1/2 + T/2t:T c ) - #(1/2), 
where ty(x) is the digamma function and T c q is the criti- 
cal temperature in the absence of impurities^. The criti- 
cal value of the scattering rate T defined by T c (r crlt ) = 
is given by r crit /T c0 = 77/27 w 1.12 within AG theory. 
However in several experiments on FeSC, e.g. Zn sub- 
stitution or proton irradiation, it is found^— that the T c 
suppression is much less than expected in the framework 
of AG theory. It has therefore been suggested that the 
s± state is not realized at all in these systems, and that 



a more conventional two-band order parameter without 
sign change (s++) is the more likely ground stat o 9 ' 10 . 

The disorder problem in these systems is substantially 
more complicated than this simple argument suggests, 
however. Even within the assumption of isotropic gaps 
on two different Fermi pockets and nonmagnetic scatter- 
ing, a much slower pair-breaking rate can be achieved 
by assuming that the scattering is primarily intraband 
rather than interband. In the pure intraband scatter- 
ing limit, Anderson's theorem applies, the system is in- 
sensitive to the sign of the order parameter, and no T c 
suppression occurs. The rate of T c suppression therefore 
apparently depends on the interplay of both intraband 
and interband scattering rates, and drawing conclusions 
regarding the superconducting state based on systematic 
disorder studies is fundamentally more difficult than in 
one-band systems. One approach to this problem has 
been to try to determine the intraband and interband 
scattering potentials microscopically for each type of im- 
purity and hostii~— , but the quantitative applicability 
of band theory to such questions is unclear. 

Here we consider the critical temperature of an 
isotropic s± two-band superconductor within the usual 
self-consistent T-matrix approximation for impurity scat- 
tering^. We perform the study both analytically in 
the weak-coupling regime and numerically in the strong- 
coupling Eliashberg framework. We find that the de- 
pendence of T c on impurity concentration is given by a 
universal form independent of impurity potentials, with 
respect to a generalized pair-breaking parameter. The 
form depends, however, on the ratio of interband to in- 
traband pairing matrix elements. Depending on the aver- 
age values of these matrix elements, we find there arc two 
possible types of s± superconductivity. The first is the 
one which has been largely discussed so far in the litera- 
ture, for which T c is suppressed as disorder is increased, 
until it vanishes at a critical value of the scattering rate. 
There is however also a second type of s± state, one for 
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which T c tends to a finite value as disorder increases; at 
the same time the gap function acquires a uniform sign, 
i.e., undergoes a transition from s± to s++. 

Model. We consider the Eliashberg equations^ for a 
two-band superconductor with a 4 x 4 matrix quasiclas- 
sical Green's function in Nambu and band space, 
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where the denote Pauli matrices in Nambu space, and 
goa and g\ a are the normal and anomalous ^-integrated 
Nambu Green's functions: 
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Here, index a runs over band indices a and b, N a ^b are the 
density of states of each band (a, b) at the Fermi level, 
and uj n = 7rT(2n + 1) is the Matsubara frequency. The 
quantities Q an and <j) an are Matsubara frequencies and 
order parameters renormalized by the self-energy S(iw n ), 
respectively, 

u an = oj n + i£oa(iu; n ) + i£o™ P ( iw n)> (3) 
4> an = Si Q (iw„) + E"" p (itJ„). (4) 

The self-energy due to the spin fluctuation interaction is 
then given by: 

£ a(iw„) = T \^<*p( n _ n ')\gop/ N P> ( 5 ) 

Si a (iw n ) = -T A Q(3 (n - n')gip/Np. (6) 

The coupling functions A a/ g(n — n') = 

/) 2 + 2 ] are expressed 
via the spectral functions (RefJ^) and constants 

\ a p- The matrix elements X a p can be positive (at- 
tractive) as well as negative (repulsive) due to the 
interplay between spin fluctuations and electron-phonon 
couplin g 15 ' 16 and strongly renormalized due to the 
nested Coulomb interaction-^. 

We use the T-matrix approximation to calculate the 
average impurity self-energy X mi P; 
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where U = U ® T3 and rii m p is impurity concentration. 
For simplicity intraband and interband parts of the po- 
tential arc set equal to v and u, respectively, such that 
(U) a/ g = (v — u)8 a p + u. This completes the specifica- 
tion of the equations which determine the quasiclassical 
Green's functions. 

Note that we have neglected possible anisotropy in 



each order parameter 



»a(b)r 



these effects can lead to 



changes in the response of the two-band s± system to 
disorder and have been examined, e.g. in Rcf. [ 



Critical temperature. T c is found by solving the lin- 
earized Eliashberg equations for the renormalization fac- 
tors Z an = u) an /u n and gap functions A Q „ = (j) an / Z a J^: 

Zan = 1 + / ] faff/l^nl 

n 

+ ttT c ^ \^ap{n - n')|sgn(w n /) /to n , (8) 

Zan^-an = / ] F a pAp n / \cj n \ 

+ ttT c ^ X a p(n-n r )Ap n >/\u n >\, (9) 
<»> n t,p 

where T a p are impurity scattering rates. 

If one inserts Eq. ((5]) into Eq. ^ and gets a set of 
equations for A Q „, it is easy to show that the impurity 
intraband scattering terms oc r oo and Tbb drop out^, in 
agreement with Anderson's theorem. From Eq. (J7J one 
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where a = (ir 2 N a N b u 2 ) / '(1 + n 2 N a N b u 2 ) and r o ( 6) = 
riimp'KNb{a)U 2 {\ — cr) are generalized cross-section and 
normal state scattering rate parameters, respectively. 
The parameter r\ controls the ratio of intra-band and 
inter-band scattering as v 2 = u 2 f]. In the Born (weak 
scattering) limit, a — > 0, while for a — > 1 the unitary limit 
(strong scattering) is achieved. From (|S4|) . we therefore 
recover explicitly the well-known but counterintuitive re- 
sult that in the unitary limit nonmagnetic impurities do 
not affect T c in an s± stat e 19 ' 20 . 

The linearized Eliashberg equations ((5J and © are 
now solved numerically, varying T and finding T c as the 
highest temperature where a nontrivial solution appears. 
Results for T c as a function of T a b, are shown in Fig lSH 
in which situation all cases with various values of a and 
rj fall on the same universal T c curve for each average 
(A) = [(A„„ + \ ab )N a /N + (X ba + \bb)N b /N] with N = 
N a + Nb. It is clearly seen that depending on the sign of 
(A) , one gets two types of T c behavior versus T a b in the s± 
scenario. For type (i), the critical temperature vanishes 
at a finite impurity scattering rate r™* for (A) < 0. For 
type (ii), (A) > 0, the critical temperature remains finite 
at T a — > 00. In the marginal case of (A) =0wc find that 
r crlt — > 00 but with exponentially small T c . Therefore, 
we have found universal behavior of T c controlled by a 
single parameter (A). 

Weak- coupling limit. To understand the origin of the 
two types of limiting behavior of T c in an s± scenario, we 
now consider the weak coupling limit assuming A Q( g (n — 
n') = A Q( 30(cjo — |^n|)0(wo~ l w nl)- In this approximation 
the calculation can be performed analytically. 
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FIG. 1: (color online). Critical temperature for various a 
and r\ as a function of the effective interband scattering rate 
r a (, for the same parameters. Note that curves for differ- 
ent sets of a and r\ overlap and fall onto one of the three 
universal curves depending on the (A). Nb/N a = 2, cou- 
pling constants for illustrative purpose are chosen for (A) > 
as (Acm, \ a b, Abo, Abb) = (3,-0.2,-0.1,0.5), for (A) = as 
(2, -2, -1, 1) and for (A) < as (1, -2, -1, 1). 
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FIG. 2: (color online). Critical temperature T c and A an (in 
units of Tco) for first, second, and third Matsubara frequencies 
n at T = 0.04T c0 , a = 0.5, t] = 1, and N b /N a = 2 and (A) > 0. 
The coupling constant are chosen as in Fig. IS1I 

We introduce the parameter 

A Q = e(co-K|)VA Q/3 7rT V ^f, (11) 

which plays the role of the pair potential in the clean 
limit. Substituting A an from Eqs. (J8|) and ([9]) and recall- 
ing r Q h/T i, a = Nb/N ai we get for |w„| < uj an equation 
for A Q similar to AG: 

A tt = A Q (A)(/ 1 -/ 2 ) + / 2 ^A Q/J A^, (12) 



where h = ^TJ2\^ n \ <UJ0 l/|w„| « ln27u;o/(7rT c ), 7 2 = 

71-71 E| W „|< W0 V(l w n| + fab + ha), A Q = E/3 and 

(A) = (A a N a /N + A b N b /N). 



In the clean limit, I 2 = h, Eq. (fT2]) reduces to A a = 
Ii^2p XapAp. Diagonalization of this equation results 
in the equation for the critical temperature I\ = 1/A , 
where A = (A aa +A bh )/2+^(A aa - A bb ) 2 /4 + A ab A ba > 
is the highest positive eigenvalue of the matrix A Q ^ . The 
critical temperature is then T c q ~ ^-^o cxp( — 1/Ao). A 
similar expression was found in&2i. The relative sign 
of the pair potential of the bands is determined by the 
off-diagonal interaction matrix elements: sgn(A a /A ) = 
sgn(A ab ). 

When A a — A b 0, nonmagnetic impurities suppress 
the critical temperature^. The critical value of the im- 
purity scattering rate for type (i) systems is given by 
ln[w /(r a b + r oa ) crit ] = (A)/(A aa A bb — A ab A ba ). 

We now focus on the case of type (ii) systems, (A) > 
0, which, to the best of our knowledge, have not been 
discussed extensively in the literature. Multiplying both 
sides of Eq. (fl"2"|) with N a , followed by a summation and 
using in the dirty limit 7 2 — ► 0, one obtains 1 = Ii(A) 
and consequently T c = ^-ojq exp(— 1/(A)). Analysis of 
Eqs.® and (j9|) shows that, for small n in the clean limit, 
A an ~ Aa, while in the dirty limit both A an and A bra 
converge to the same value, A an — > Ar^oo, that is the 
s ++ state is realized. If the initial state corresponds to 
s±, a transition s± —> s++ at a finite concentration of 
impurities must exist. 

There is a simple physical argument behind the s± — > 
s ++ transition. With increasing inter-band disorder, the 
gap functions on the different Fermi surfaces tend to the 
same value. A similar effect has been found in Refs.— ^ 
for a two-band systems with s ++ symmetry, and in Ref.— 
discussing node lifting on the electron pockets for the 
extended s-wave state in FeSC. 

To demonstrate the transition explicitly, we calculate 
A an for n = 0, 1, 2 at T = 0.04T c o and show the results 
in Fig. [2]for a particular choice of \ a p with (A) > 0. For 
this parameter set, T c q « 40K. Both order parameters 
A a (h)„ converge to Ap-s-oo for large disorder, while the T c 
suppression quickly saturates. The transition s± — > s ++ 
provides a possible explanation for the observed much 
weaker reduction of the critical temperature than the 
naive application of the AG formula. 

Another important consequence of the transition s± — > 
s_l__l_, relevant to experiments in pnictides, is gapless su- 
perconductivity as one of the gaps vanishes. The den- 
sity of states Ntoti 1 -^) = — J2 a I m 5oa( w )/ 7r is shown in 
Fig. [3l[a) for a type (ii) case. With increasing impurity 
scattering rate, the lower gap is seen to close, leading to 
a finite residual A tot (w = 0), followed by a reopening of 
the gap. A similar behavior is reflected in the temper- 
ature dependence of the penetration depth, [Fig. EJb)], 
which varies in the clean limit with activated behavior 
controlled by the smaller gap, crossing over to T 2 in the 
gapless regime, to a new activated behavior in the s++ 
state in the dirty limit. Figure EJb) should be compared 
to similar works, where the effect of scattering on the 
T-dcpcndent supcrfluid density was calculated for a two- 
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(a) 




FIG. 3: (color online), (a) Density of states Ntot(u)/N vs. 
r a /T c o and ui/T c a for (A) > and impurity parameters a — 
0.5, T) = 1, N b /Na = 2, N = N a + N b . (b) Total superfluid 
density 1/(lj p \l) 2 vs. F a /T c o and T/T c o where lo p is the total 
plasma frequency and Al is the London penetration depth. 



band s + + statei^^^. 



Conclusions. We have shown that in two-band models 
with an s± ground state, T c has a universal dependence 
on the impurity scattering rate which can be calculated 
explicitly in terms of the interband to intraband impurity 
scattering rate ratio. We demonstrated that s± super- 
conductivity may be quite robust against nonmagnetic 
impurities, depending on the ratio of interband to intra- 
band pairing coupling constants, and may even display 
a transition to an s ++ gap structure with increasing dis- 
order, which will manifest itself in thermodynamic and 
transport properties. 
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I. SUPPLEMENTARY ONLINE MATERIAL FOR THE ARTICLE "DISORDER INDUCED 
TRANSITION BETWEEN s± AND s++ STATES IN TWO-BAND SUPERCONDUCTORS" 

Multiband system in the weak coupling approximation. 

The Eliashberg equations in the general form on the imaginary Matsubara axis are 



= *t£E x ai (n - n>) m _ + £ f 



W«n — U) - 



J) 2 

'in 1 T%n 



We consider the weak coupling limit (the so-called 90-modcl), which corresponds to: A a p{n — n') = \ a p<d(oJo 
|w n |)0(wo — l w nl)- in this case, the term with X a i{n — n') vanishes, 

4>an = 6(^0 - |w„|)^A M 7rT ^ = z h ^ f . 



T f„' + #n' 1 V ^in + $n 
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Zan = 1 + £ 
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Note that intraband nonmagnetic impurities scattering rate T aa drop out from this equation according to Anderson's 
theorem. Let us consider for simplicity T = T c . We would also like to split the gap functions into two parts: part A 
undergoes a strong impurity scattering between different parts of the Fermi surface separated by large wave vectors 
(say, electron and hole bands) with oc UiNj , while the part B undergoes weak scattering with 7^ oc U2Nj. In 
this case, Eqs. (|S1|) for a € A become 

A a „ = 0(cj o - |w„|) V" A ai 7rT c V* 

I \UJ n /\< : UJO 

iY^f A jn - A an A/ n - A an 

+ £ r ^ | Wn | + E ^ | | - (S2) 

jeA 1 (es 1 

For \u> n \ < wo, the solution of this equation can be written as 

r _ J2jeA r aj^jn + E/es7a/A in +/ |w n | 



where 

/ a =7rT c VA ol V (JVi ~ ^ /2^T c0 » 1). 



(S3) 
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Since r a j oc 1/4 TVj, for u\ y oc from Eq. (|S3|) wg see that the gap function coincides with A an oc Yljei Njlj- On 
the other hand, the gap functions belonging to the other part (weakly coupled bands) can have different signs. 
T c dependence on the effective interband scattering rate. 
One finds that f ab(ba) = T a ( b )f(a,r)) 
with 

f{ *> = - ri ^ (w , r —&> (S4) 

where a = (it 2 N a N b u 2 ) / (1 + it 2 N a N b u 2 ) and r a ( b ) = ni mp 7riV b ( a )ii 2 (I — cr) are generalized cross-section and normal 
state scattering rate parameters, respectively, and ni mp is the impurity concentration. The parameter 77 controls the 
ratio of intraband and interband scattering as v = un. In the Born limit, a — > 0, while for a — > 1 the unitary limit 
is achieved. Function /(a, 77) is limited from above by 1, therefore, in the unitary limit nonmagnetic impurities do 
not affect T c . From Eq. (|S4[) . we therefore recover explicitly the well-known but counterintuitive result that in the 
unitary limit nonmagnetic impurities do not affect T c in the s± state. 

Results for T c as a function of the pairbreaking parameter F a (proportional to impurity concentration ni mp ) are 
shown in Fig. ISlTa). For illustrative purposes, the coupling constants A Q( g are chosen the same as in Fig. 1 of the 
main text. Note that the strongest suppression is generally found for pure uniform scattering, 77 = 0, in the Born 
limit (a 0) 

and that in the opposite limit of pure intraband scattering, u = 0, we have rj — > 00, so that there is no pairbreaking 
since r ab — » 0. The similar situation takes place for the strong unitary limit. It is clearly seen that depending on the 
sign of the average (A) = [(A 00 + X a b)N a /N + (\ ba + \ bb )N b /N] with N = N a + N b , one gets two types of T c vs. r a 
behavior in the s± scenario. For type (i), the critical temperature vanishes at a finite impurity scattering rate F^ rlt 
for (A) < 0. For type (ii), the critical temperature remains finite at T a — > 00. In the marginal case of (A) = we find 
that T™* — > 00, but with exponentially small T c . Fig. ISlfb) is equivalent to Fig. I of the main text and shown here 
for easier comparison with the panel (a) . 
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FIG. SI: (color online). Critical temperature for various a and rj as a function of (a) the impurity scattering rate r a and (b) 
the effective interband scattering rate T a b for the same parameters. Nb/N a = 2, coupling constants are the same as in Fig.l of 
the main text. 



